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> 1 ABSTRACT 

i> : 

m . 

Neutron production via 4 He breakup and p(p, nix + )p is considered in the 
innermost region of an accretion disk surrounding a Kerr Black Hole. These 
reactions occur in a plasma in Wien equilibrium, where (radiatively produced) 
pair production equals annihilation. Cooling of the disk is assumed to be due 
to unsaturated inverse Comptonization of external soft photons and to the en- 
ergy needed to ignite both nuclear reactions. Assuming matter composition of 
90% Hydrogen and 10% He, it is shown that, close to the border of this region, 
neutron production is essentially from 4 He breakup. Close to the horizon, the 
.^J \ contribution from p(p, n7t + ) p to the neutron production is comparable to that 

from the breakup. It is shown that the viscosity generated by the collisions of 
the accreting matter with the neutrons, may drive stationary accretion, for ac- 
cretion rates below a critical value. In this case, solution to the disk equations is 
double-valued and for both solutions protons overnumber the pairs. It is claimed 
that these solutions may mimic the states of high and low luminosity observed in 
Cygnus X-l and related sources. This would be explained either by the coupling 
of thermal instability to the peculiar behavior of the viscosity parameter a with 
the ion temperature that may intermittently switch accretion off or by the im- 
possibility of a perfect tuning for both thermal and pair equilibrium in the disk, 
a fact that forces the system to undergo a kind of limit cycle behavior around 
the upper solution. 



Subject headings: accretion disks-nuclear reactions-radiative transfer-pair pro- 
duction 
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1. Introduction 



Accretion disks are presently thought as the scenario for hydrodynamical flows in close 
binary systems, galactic nuclei and quasars. Despite this fast growing recognition, models 
which allow for some particle processes, taking into account the detailed nature of the flow, 
as well as realistic processes that may act as the source for the required viscosity, are still 
lacking. Concerning particles processes, most of the theoretical work on the field has focused 
on the production of positron-electron pairs. Steady state scenarios, under the assumption 
of production-annihilation equilibrium, have been first tackled by Bisnovatyi-Kogan et al. 
(1971), who found that, if particle-particle dominates pair creation, equilibrium is only pos- 
sible for lepton temperature below 20mc 2 , Pozdnyakov et al. (1977), Stoeger (1977), and by 
Liang (1979), who considered pair production dominated by gamma-gamma interactions, in 
a plasma under Wien equilibrium, finding multivalued solutions to the disk structure. These 
results, however, are applicable only to a small region of parameter space or to a subset of 
important reactions. More general studies have been subsequently carried out by many au- 
thors (Lightman 1982; Svensson 1982, 1984; Zdziarski 1984; Takahara and Kusunose 1985; 
Lightman and Zdziarski 1987; White and Lightman 1989), which contributed to a better 
understanding of radiative processes in very hot astrophysical plasmas, opening the way to 
the understanding of the topology of the disk equations solution, which at fixed r in the M £ 
plane forms an S-shape kink corresponding to its multiple valued nature (Kusunose 1990; 
Liang 1991; Bjornsson and Svensson 1991; Bjornsson and Svensson 1992; Minishige 1993; 
Kusunose and Minishige 1992, 1995, 1996). This multivalueness consists of three branches, 
two of them being hot (one pair dominated and the other pair deficient) and the third being 
cool and pair deficient. The hot pair dominated branch is very promising to explain some 
especial features in some black hole candidates, like the MeV bump of Cygnus X-l, the 
radio plasmoid bipolar outflow of the jet sources, and the alleged annihilation line features 
(Ling et al. 1987; Liang 1998; Mirabel et al. 1992; Mirabel and Rodriguez 1994; Hjellming 
and Rupen 1995; Paul et al. 1991). It may happen that, under the conditions prevailing 
in the innermost parts of some systems, the time needed for the matter to cross these re- 
gions is comparable to, or even less than, the time the electrons need to extract energy from 
the protons, through collisional energy exchange. These regions will grow in the direction 
perpendicular to the plane of the disk, becoming geometrically thick (swelling of the disk 
may also occur for systems with Eddington or super Eddington luminosities (Abramowicz 
et al. 1980; Paczynski and Abramowicz 1982), not considered in this paper). From angular 
momentum transport point of view, collective processes are efficiently operative to generate 
viscosity (due to the protons), but with a deficient cooling (due to the electrons). This drives 
protons and electrons out of thermal equilibrium, the ion temperature being much greater 
and close to the virial one. For a disk around a maximally (synchronously) rotating Kerr 
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black hole the inner radius of the disk may reach the horizon, with the ion temperature 
approaching the mass of the proton. A two temperature geometrically thicker and optically 
thin disk model was first developed by Shapiro et al. (1976) to explain the hard X-ray from 
Cygnus X-l. These authors, however, have not considered radial advection of energy and 
entropy which, in that situation, are very important and the local approximation breaks 
down. Advection of energy in optically thin disks was first considered by (Ichimaru 1977), 
in the context of the bimodal behavior of the X-ray spectra from Cygnus X-l, (Liang and 
Thompson 1980) who considered advection in an optically thin, two-temperature disk, and 
by (White and Lightman 1990) who studied advection of energy and pairs in optically thin 
disks. (Paczynski and Bisnovatyi-Kogan 1981; Muchotrzeb and Paczynski 1982; Abramow- 
icz et al. 1988; Honma et al. 1991; Wallinder 1991; Chen and Taam 1993), have unraveled 
the basic physics of advection in accretion flows with their studies on the properties and 
structure of optically thick disks. However, only recently, advective cooling in optically thin 
accretion flows has been recognized thanks to a systematic work by Narayan and Yi (1994, 
1995), Narayan et al. (1995a,b), Abramowicz et al. (1995), Narayan et al. (1995a,b), Narayan 
et al. (1996), Lasota et al. (1996a,b), and Chen et al. (1995). Neglecting pairs, these au- 
thors have studied the topology of the solution in the ME plane, showing the existence 
of a maximum accretion rate, above which no steady state solution is allowed, and that, 
below it, there are two optically thin solutions, namely the radiative cooling dominated and 
the advective cooling dominated states. A general description of flows in accretion disks, 
taking advection into account and neglecting pairs, was given by Chen et al. (1995). They 
have shown that, at a given radius, there exist exactly four physically distinct types of 
accretion disks. Two of these correspond to values of a, the viscosity parameter, smaller 
than a critical value, and the other two, to values of a greater than this critical value. In- 
clusion of pairs have been considered by Bjornsson et al. (1996), Kusunose and Minishige 
(1996), Esin et al. (1996, 1997), who have shown that pairs do not modify the topology of 
the solution, besides being negligible for a < 1. It should be remarked that, though the 
advective model constitutes an improvement as far as the standard Shakura and Sunyaev 
disk model is considered, till the moment, no global solution exists for the disk (Bjornsson 
et al. 1996). At this point, we must realize the importance of the viscosity parameter and, 
yet, we don't know the physical process that may be the source for such a viscosity. Since 
the seminal paper by Shakura and Sunyaev (1973), a lot of mechanisms have been proposed 
to account for viscosity: shear turbulence generated by the Keplerian rotation of the disk 
(Shakura and Sunyaev 1973; Zel'dovich 1981; Dubrulle 1990; Zahn 1991), turbulence driven 
by convection (Lin and Papaloizou 1980; Tayler 1980; Bisnovatyi-Kogan and Blinnikov 1977; 
Shakura et al. 1977; Ryu and Goodman 1992; Cabot and Pollack 1992; Stone and Balbus 
1996; Meirelles 1991a, 1993, 1997), tangled magnetic fields sheared by differential rotation 
(Lynden-Bell 1969; Shakura and Sunyaev 1973; Eardley and Lightman 1975; Ichimaru 1977; 
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Coroniti 1981), angular momentum transport by waves (Papaloizou and Pringle 1984a,b; 
Vishniac and Diamond 1989; Vishniac et al. 1990), Velikov-Chandrasekhar magnetic insta- 
bility (Balbus and Hawley 1991, 1992; Balbus et al. 1996; Vishniac and Diamond 1992), 
ion viscosity and neutron viscosity (Guessoum and Kazanas 1990; Meirelles 1993), radiative 
viscosity (Loeb and Laor 1991). None of these processes, however, is immune to criticism: 
concerning shear turbulence, disks satisfy Rayleigh's criterion for stability; convection may 
transport angular momentum inwards rather than outwards; magnetic field may be removed 
by magnetic buoyance and ohmic dissipation; concerning wave turbulence, it seems that the 
most stable mode have low wavenumbers, giving rise to structures of the order of the size 
of the system, being very sensitive to the Coriolis force which hinders the appearance of 
smaller structures; the Velikov-Chandrasekhar instability is highly dependent on the radial 
structure of the azimuthal component of the magnetic field, and can only occur if the Alfven 
velocity is of the order of the Keplerian one, or if the typical scale of the magnetic field is 
smaller than the Keplerian one; as far as ion viscosity is concerned, besides needing high 
temperatures (the system needs another process to achieve these high temperatures), any 
small magnetic field present in the flow will decrease the strength of the viscosity by orders 
of magnitude; neutron viscosity needs high temperature to ignite nuclear reactions capable 
to produce neutrons, needing, therefore, some other process to heat up particles till these 
high temperatures; Keplerian thin disks cannot be supported by radiative viscosity, because 
the required energy density of the photon gas should be one or two orders of magnitude 
larger. 

In view of all these uncertainties concerning the hydrodynamics of the flow, as well as the 
angular momentum transport in accretion disks, we would like to consider neutron viscosity, 
assuming neutron production through Helium break up and through pion production due to 
proton-proton collisions. The threshold for Helium breakup is about 20 MeV and about 290 
MeV for the production of pions (and neutrons) in the proton-proton collision. Close to the 
horizon, the energy in the rest frame of two colliding protons is about 2 GeV, which means 
that, at least theoretically, these reactions are energetically viable. The real difficulty is to 
find conditions in parameter space, if they exist, such that the plasma is no longer ruled by 
electromagnetic interactions alone and nuclear interactions start playing a role. Effectiveness 
of strong interactions depends not only on the value of ion temperature itself but also on 
the electron temperature. This is so because electromagnetic interactions are restricted to 
the Debye sphere and the number of electrons within it decreases with electron temperature. 
High electron temperature will favor, not only the production of electron-positron pairs, but 
also the possibility of nuclear reactions. One has to realize we are here dealing with a very 
intricate situation highly dependent on the viscosity. It is well known, that for a given density, 
electron and proton temperatures, drift time decreases with increasing viscosity parameter 
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while t ep , the electron-proton collision time, has inverse behavior. However, physical variables 
in the disk are very sensitive functions of the viscosity parameter and we have to make a 
self-consistent calculation to find a region in parameter space where the nuclear reactions 
effects are maximized. 

Account for nuclear reactions in the accretion disks may have interesting observational 
consequences such as the lines spectral features of Hydrogen-Helium plasmas at high tem- 
peratures and 7-ray lines production. From the theoretical point of view, Helium breakup 
and pion production reactions produce neutrons whose collisions with the accreting matter 
may be a source for the viscosity in the disk, which, in turn, will depend upon electron tem- 
perature, proton temperature and accretion rate. This may have interesting consequences 
as far as the topology of the solution is concerned. This is a problem we want to tackle in 
a future work. However, in the present paper, it should be stressed, our main concern is to 
find out if neutron production through proton-proton collisions is viable in the innermost 
region in accretion disks, and if its collisions with the accreting matter can drive accretion 
on. We will, therefore, be interested only in order of magnitude estimates. Keeping this on 
mind, we will make some reasonable approximations, mainly on the radiative transport, pair 
production and on the hydrodynamics, which, at the right moment, will be justified. 

A potentially promising application of the afore mentioned reactions is the modelling of 
some X-ray systems, like Cygnus X-l, that exhibit multimodal behavior. In a previous paper 
(Meirelles 1993), it was shown that although these reactions, in steady state accretion, can 
not supply the disk with enough neutrons to make their collisions with the accreting matter 
the main source of viscosity, they may have strong implications as far as the temporal 
behavior of the disks is concerned. We have considered, however, only the production of 
neutrons through the 4 He breakup reaction 

4 He + 28.296Me\/ -> 2p + 2n, (1) 

for ion temperatures greater than 3 MeV. Assuming that a steady accretion can be achieved 
with the drift time equal to the nuclear reaction time, we were able to show the existence of 
a critical accretion rate, below which there is no steady state accretion onto the hole. Above 
the critical accretion rate it is possible, under special circumstances, for the disk to choose 
between the two states of steady accretion. 

This kind of procedure has to be criticized on the grounds that equality between drift 
and nuclear reaction times is a very stringent constraint and, at the temperature range 
considered, electron-positron pair production should be taken into account. Besides, as we 
approach the hole, ion temperature can be much greater and even exceed the threshold 
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temperature for the reaction 

p + p + 290MeV ^ p + n + ir . (2) 

We claim that accounting of both reactions in the inner parts of an accretion disk to- 
gether with allowance for thermal instability may explain the transitions observed in Cygnus 
X-l, between the states of high and low luminosity. 



2. Observational aspects: a brief review 

There is a widely accepted suspicion that the unseen compact object in the binary sys- 
tem Cygnus X-l is a black hole (Liang 1998). Support for this suspicion is found on detailed 
analyses carried on radial velocity measurements as well as on recent analyses based on spec- 
trum synthesis (Herrero et al. 1995) both of which give a mass of about 10 M , comfortably 
larger than the 3 M upper limit of a neutron star mass. The luminosity is sub-Eddington, 
but greater than 0.01 L E dd (Liang and Nolan 1984). Cygnus X-l, being one of the most 
studied sources in the sky, especially in the hard X-rays, some of its properties (characteris- 
tic hard X-ray spectrum, episodic emergence of an ultra soft component, anticorrelated soft 
and hard X-ray transitions, chaotic variability down to milliseconds, persistent gamma-ray 
tail above an MeV, episodic gamma-ray bumps at a few hundred keV-MeV, persistent radio 
emission and radio flaring correlated with X-ray transitions and Low-Frequency Quasi Pe- 
riodic Oscillations (QPO)) have been accepted as evidence for black hole candidacy (Liang 
1998). Observations usually find the system in one of the two states: the hard state and the 
soft state (Oda et al. 1974; Liang and Nolan 1984; Tanaka and Lewin 1995). Most of the 
time, Cygnus X-l is found in the hard state, where its soft X-ray (2-10 KeV) is relatively low 
and the hard X-ray spectrum is hard. After being in this state for a few years, the system 
undergoes a transition to the soft state, increasing its soft X-ray flux by a factor of about 
10, softening its X-ray spectrum. It remains in the soft state from weeks to months before 
returning to the hard state. The transition between the two states lasts from less than a day 
to more than a week (Wen et al. 2001). In the hard state, Cygnus X-l has a power law spec- 
trum characterized by a photon index ss 1.4. During the hard state the soft X-ray spectrum 
below 10 KeV is often a simple continuation of the hard X-ray power law, with increasing 
flattening below 3 KeV. In that case, the entire X-ray continuum is likely produced by a 
single hot component. However, during some episodes of Cygnus X-l hard state, the best 
fit model still requires a small black-body component of temperature a few tenths of KeV 
on the top of the power law, presuming the existence of another region of lesser temperature 
(Liang 1998). This is highly suggestive of a soft photon source, as expected in the inverse 
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Compton model. In addition there is often evidence of the FeK-edge absorption at w 7 KeV 
(Liang 1998; Esin et al. 1998). When Cygnus X-l is in the soft state, its spectrum switches 
to one that is dominated by an ultrasoft component, with a temperature ~ 0.3 — 0.4 KeV. 
The spectrum above 10 KeV becomes much softer than in the hard state, with a variable 
photon index « 1.9 - 2.5 (Gierlinski et al. 1997; Cui et al. 1997; Dotani et al. 1997). It 
exhibits, also, pair annihilation features in the region 500 < E < 1000 KeV (Liang 1993a). 



3. The model 

The region we shall be concerned extends away from Ri, inner radius of the disk, to 
approximately 10 to 15 Ri, where R^ -5Rs, Rs being the gravitational radius. 

We shall assume conditions such that the ion temperature, T i: is much greater than T e , 
the electron temperature, with pressure given by the ions. In the outermost parts of the 
disk, matter composition is Hydrogen (nine tenth) and Helium (one tenth). In its way down 
to the hole, He depletion starts and the neutrons produced gradually contributes to the 
viscosity. As the incoming matter comes closer to the hole, charge exchange reactions start 
to contribute to the neutron production. Concomitant to these reactions, electron-positron 
pairs are radioactively produced. We shall assume production-annihilation equilibrium. We 
shall also assume an external source of soft photons that are continuously impinging upon the 
disk. Part of these soft photons succeed being upscattered in energy, reaching the Wien peak 
where they interact and produce pairs. The remaining are upscattered, but leave the disk 
before reaching the Wien peak. Both nuclear reactions we are considering are endothermic, so 
they lower the plasma temperature by taking the energy they need to ignite. For the Helium 
breakup, the number of the resulting particles is greater than the number of the interacting 
particles. These particles have to be thermalized and so act as an additional source of 
cooling. The pion produced by the charge exchange reaction rapidly decay, producing a 
positron that has to be thermalized and a neutrino that leaves the system, cooling the disk 
even more. Our model has axial symmetry, the azimuthal velocity being keplerian, and 
hydrostatic equilibrium holds in the direction perpendicular to the plane of the disk, with 
the pressure given by the nucleons. The disk is heated by viscous processes and locally cooled 
by radiation and by nuclear reactions. This innermost region is surrounded by a region in 
which matter is completely depleted of He, due to breakup. Out of this neighboring region, 
we assume the existence of some unknown process that heats the gas and makes He breakup 
viable. Since the nucleon temperature should be close to the virial one for charge exchange 
reactions to occur, the disk will be geometrically thick, and advection should be important. 
However, the more important advection is, the less important is the radiative cooling, and 
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the greater the nucleoli temperature will be, a fact that favors the onset of nuclear reactions. 
Contradictory as it may seem, we don't know by sure the effect on the nuclear reaction rate 
by the neglecting of advection. Besides this, Esin et al. (1997) have studied the behavior 
of advection as a function of the accretion rate and of the viscosity parameter, concluding 
for its importance of for small values of m (in Eddington units). As rh increases, from 0.01 
to a critical value that depends on the viscosity parameter, the advection zone becomes 
very luminous. The advection zone shrinks in size with further increase of the accretion 
rate, and for sufficiently large values it disappears, with the thin disk extending down to 
the marginally stable orbit. In this paper we are interested in order of magnitude estimates, 
and accounting of advection would make our calculations very complex. Meirelles (1991a) 
has taken conductive transport into account in the inner geometrically thick region of a two 
temperature soft Comptonized accretion disk in Cygnus X-l, showing that this process is 
much more important than radiative cooling. For a matter of consistency, were advection 
included, we should also have included conduction, which would render the problem much 
more complicated. 



4. Disk Equations 

Hydrostatic equilibrium in the vertical direction reads 

P=^Q 2 H 2 , (3) 
3 

where p is the density, Q is the Keplerian angular velocity and H is the semi-scale height of 
the disk. P, the total pressure, is mainly given by the pressure due to the protons, i.e., 

P = N 10 11 kTi , (4) 

N being the proton number density, k the Boltzmann constant. Ion temperature is in units 
of 10 11 and electron temperature in units of 10 9 . Unless otherwise stated, physical variables 
are expressed in the cgs system of units. Combining eq. (3) and (4) we obtain 

H = 4.93 x 10 9 IT 1 Ti - 5 . (5) 

Mass conservation, together with conservation of angular momentum and the definition 
of the viscosity parameter yield 

^{l = aPH, (6) 

4 7T 

which implies for the column density 

E = 1.96 x 10~ 20 a' 1 T' 1 Q M S , (7) 
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a being the viscosity parameter, M is the accretion and S is a function that depends on 
the boundary condition imposed on the torque at the inner radius. For Q + , the heat flux, 
viscously generated, going into the disk at a distance R from the central object, 

Q + = 3/8vr MVl 2 . (8) 

In the two-temperature regime, energy is collisionally transferred from protons to elec- 
trons and pairs at a rate, corrected for the inclusion of neutrons, given by (Guilbert and 
Stepney 1985; Spitzer 1962) 

Q~ = 1.72 x 10 28 p 2 (Tj - T e ) (1 + 0.41 T e a5 ) (1 + 2 z + y + ) (1 + y n Y 2 T~ z ' 2 H , (9) 

where y n is the total neutron abundance, y + is the abundance due to pion production and 
z is the positron density to proton density ratio. The electrons and pairs, in turn, loose 
energy through unsaturated inverse Comptonization of externally supplied soft photons. 
Like Shapiro et al. (1976), we shall assume these photons are copiously supplied and that 
the Kompaneetz y parameter is constant and equal to 1 (Rybicki and Lightman 1979), i.e., 

y = (0.674T e ) 1 rj = 1, (10) 

where % — 1 for T e < 5.9 and 2 otherwise, j = 1 for r es < 1 and 2 otherwise. We shall 
assume that the intensity of the soft photons coming from the external source overwhelms 
the intensity due to photons internally produced (bremsstrahlung) . With the inclusion of 
neutrons, the electron scattering depth changes to 

T es = ^Y,{l + 2z + y + ){l+y n y 1 , (11) 

a T h being the Thomson cross section for electron scattering and m p the proton mass. Clearly, 
this way of treating the radiative problem in the disk relies heavily on the existence of 
an external soft photon source, which we assume, from the very beginning, as granted. 
We should emphasize, since our main concern is not the radiative problem, we adopt this 
procedure for practical reasons: the spectrum of these soft photons, after being upscattered 
in energy, through unsaturated inverse comptonization in the inner region, reproduces fairly 
well the observed spectrum of Cygnus X-l, in the 8-500 KeV region (Shapiro et al. 1976). 

For the moment we postpone the discussion about the cooling of the disk due to breakup 
and pion production. We next discuss pair production in the disk. 

5. Pair Production 

As we have emphasized before, our main concern in this paper is find out the possibility 
of occurrence of nuclear reactions in the inner region of the disk. We, therefore, construct a 



-10- 



reasonable and realistic environment of photons and pairs, in which the nuclear reactions take 
place. As we have done in the case of radiative transport, we shall make some simplifications 
when treating the pairs. We assume pairs are created mainly by Wien photons interacting 
with Wien photons. However, instead of writing down the formal pair equilibrium equation, 
we shall adopt a procedure similar to Liang (1979) paper. This essentially consists in relating 
the photon chemical potential to the soft photon flux. This kind of procedure is valid as 
long as the plasma is in Wien equilibrium (Svensson 1984), which implies r 77 > 1. Below, 
we briefly review this procedure. 

For kT e < m e c 2 , f u , the photon occupation number, obeys the Kompaneetz equation. 
In that limit 

/„ = 24 x~ s (l+x + x 2 /2\ + x 3 /3! + x 4 /4\) , (12) 

where \x v is the photon chemical potential and x — hvj kT e . The occupation number / + _ 
for positrons and electrons may be written, respectively, 

/+ = e "+e- E +, (13a) 

f_ = e»-e- E -. (136) 
Since the number of particles is conserved in the reaction 

7 + 7 -» e + + e~ , (14) 

the chemical potential should satisfy 

2 /i 7 = /!+ + /i_ . 15 

For Cygnus X-l and related sources n ~ 1, which implies s — 4. The radiative cooling will 
be 



7T / — —dv. (16) 

Jv c 



Now, using eqs. (12), (13), (14), (15) and (16) we obtain 

' Jc 2 F r 
\kT e f 276 7T 



\ ! ( „ii'L ) • (17) 



where J is given by 



./ \-r,n,-kT. ^i(^y-) , (18) 



K 2 being the modified Bessel function of 2nd kind and we have substituted in the slowly 
varying function In (k T e /h v ) for canonical values appropriate for Cygnus X-l. 
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lt is worth mentioning that the radiative cooling now will be given by 

F r = Q + -Q n , (19) 

where Q n is the nuclear cooling. Adopting the procedure as we did to treat radiative cooling 
and pair production we have overestimated both processes. Since the energy production is 
constant, this implies in underestimation of the nuclear reaction rate. 



6. 4 He, pion and neutron production 

In the chain of reactions leading to the 4 He breakup 

18.35 MeV + p + 4 He — > 3 He + D 

20.578 MeV + p + 4 He -> 3 He + p + n 

19.844 MeV + p + 4 He -> 3 H + p + p 

5.494 MeV +p + 3 He -> D+p + p 

0.739 MeV + p + 3 H ^ 3 He + n 

2.224 MeV + p + + p + n 

the one that takes longer to occur is the first one. Once this reaction has ignited, the other 
follow very rapidly. So, we assume the breakup dominated by this reaction (Guessoum 
and Dermer 1988) and use a prescription (Gould 1982, 1986; Guessoum and Gould 1989; 
Guessoum and Kazanas 1990) to write the neutron reaction rate for this reaction, 

R b = 5.67 x 10~ 16 T~ - 5 e" 2 - 56/Tl N 2 . (20) 

This implies a neutron abundance, due to this reaction, given by, 

y b = imxl0 9 T-°- 5 e- 2 - 56 / T ^ 2 R 2 H- 1 (1 + y n y l . (21) 

The total neutron abundance y n will be given by 

y n = y+ + y b . (22) 

Using eq. (20) we can express the nuclear cooling due to this reaction as 

Q ei = 2.84 x 10 20 (1 + 1.37 TO T~ l S 2 (1 + y n y 2 . (23) 



- 12 - 



For the reaction 

p + p^p + n + n + 

we shall use the results by Engel et al. (1996), who have carried out a very detailed study 
of the p (p, nn + )p and p (p, pn°)p reactions using a fully relativistic Feynman diagram tech- 
nique. Their calculations have been carried out under the prescription of the one boson 
exchange (OBE) model, with allowance for the inclusion of both nucleon and delta isobar 
excitations in the intermediate states, as well as for the exchange of n, a, p, and u mesons. 
Accounting for the exchange of a and oo makes their results very reliable close to the thresh- 
old, where these contributions presumedly dominate. Most of the parameters of the OBE 
model are determined by fitting to the N-N scattering data over the energy range of 300 
MeV to 2 GeV. A fitting to their data for the total cross section leads to 



o vv = 10 f(E) . (24) 



where 



0.0371785E - 14.781764 , 290.0 sC E < 378.8Me1/ 

/ (E) = { -4.96247 + 0.0148283E - 0.00000839163E 2 , 378.8 «C E ^ 1060Me\/ 
-0.000064893617E + 1.430787 , 1060 ^ E ^ 2000Me\/, 



where E is expressed in MeV and the cross section in mb. The reaction rate for this reaction 
will be written as 

R+ = 0.5N 2 i — - / a pp e- w wdw, (25) 

\ m p / Jw s 

w is the energy in units of kTi and w s is the threshold energy for this reaction. Then, y + , 
the contribution of this reaction to the neutron abundance, follows, straightforwardly, 

y + = 9.24 x 10~ 4 R 2 S 2 (^) ° ' t(w) M^ 1 H' 1 (1 + y n y l , (26) 

V m p J 

where t(w) is the integral defined in eq. (25), now in units of 10~ 27 mb. Figure 1 shows t (w) 
as a function of T im 

To write down the cooling of the disk due to this reaction, we remind ourselves that the 
pion has a mean life of about 2.6 x 10~ 8 s, decaying through 

tt + -> fi + + + 34 MeV . 

The in turn, has a mean life of about 2.2 x 10~ 6 s, decaying through 

/x + -> u e + e + + + 105 MeV , 
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which allows us to write for the cooling 

Q pp = 1.23 x 10 29 (2.40 x 1(T 2 + 2.07 x 10~ 7 T e ) £ 2 Tj ' 5 1 (w) H' 1 (1 + y n y 2 . (27) 

Finally, we can verify if neutrons can respond for the viscosity in the disk. In order to 
do so, we write for the neutron viscosity (Weaver 1976) 

7] = l/3p n v n \ n , (28) 

where p n , v n , and A n are respectively neutron density, neutron velocity and neutron mean 
free path. Taking the average for a Maxwellian distribution, we obtain (Bond et al. 1965) 

/ k T \ 5 

^ = 1.08j/ n — M e n , (29) 
\mn J 

£ n being the effective mean free path given by 

in = , (30) 

with A n given by 

A. = ^— ■ (31) 

a n p{l + y n ) 

The cross section for neutron scattering, averaged over a Maxwellian distribution, is (Gammel 
1963) 

a n = 5A7T-°- 85 b . (32) 
Finally, using the definition of the viscosity parameter a together with eq. (29) we arrive at 

a 2 = 3.65 x 10~ 2 M l7 y n Q T' 185 (1 + y n y l - 10~ 5 M 2 7 fi 2 T^ 3 ' 7 , (33) 

the accretion rate expressed in units of 10 17 gs _1 . 

We are now in a position to ask ourselves if neutrons collisions with the accreting 
matter can indeed act as a source for the viscosity in the innermost regions of an accretion 
disk. We shall see that previous results (Meirelles 1993) relied heavily on the constraint 
imposed assuming equality between nuclear reaction time and dynamical time. Dropping 
this assumption leaves practically no restriction as far as the accretion rate is concerned. To 
not be burdened by additional complexities let us unravel things a little bit by neglecting, for 
the moment, the contribution to the neutrons due to pion production. Under this procedure, 
combination of eqs. (5), (7), (21) and (33) yields 

y n 2 (3.65 x 10~ 2 MnilXi 1 - 15 - 10~ 5 M 2 7 fi 2 Xp a7 j 



-14- 



-y n 10~ 5 M 2 7 Q 2 Tr - 7 - 0.458 e" 2 - 56/Ti = . (34) 
Now specializing for r = 10, we obtain 

Mi7< 5.77 Ti 1 - 85 . (35) 

This, however, is not the final result since we have to investigate the effects of cooling and 
pair production. Figure 2 shows the neutron abundance y n as a function of T, for r = 10 
and M 17 = 1.0. 

In Fig. 3 we present the viscosity parameter a as a function of T, in the same conditions 
of Fig. 2. We see that the viscosity parameter may reach very high values. The maximum 
a occurs for % ps 3.0. 

We finally reduce our system of equations for the disk to only two equations involving 
T e and T;. The first of them is the thermal equilibrium equation and reads 

(l + 1.37T l )e" 2 - 56 ^a- 2 T r 3 (l + y„)- 2 

+ 12.14 (1 + 0.41 T e a5 ) a' 1 T> 5 T e 2 ' 5 (1 + y n )~ l - 10.92 = , (36) 
and the second is the pair equilibrium equation, 

A 2 - 1 - 252.5 T 2 (l + 0.41T e a5 ) 2 e~ n - 87/Te T e - 10 = , (37a) 

where A is given by 

A = 6.26aT J T e - 1 (l + y n ) 1 - 2/n , (376). 

a and y n are given respectively by eqs. (33) and (34). 

A glance at Fig. 4 reveals the existence of a solution for the two-temperature soft photon 
Comptonized accretion disk with pairs and viscosity generated by neutron collisions with the 
accreting matter, neutrons supplied by 4 He breakup. 

To check for the consistency of our solution we also plot (1 + 2z) in Fig. 5, both in the 
thermal equilibrium and pair equilibrium situations. 

7. Highlighting pion production 

We have seen in the previous section, neglecting contribution from pion production 
to the neutrons in the disk, the existence of a solution to the disk equations at r = 10, 
Tj ps 2.5. Were pion production included, the results would not change so much, since at 
this temperature the contribution would be negligible. However, at that temperature, the 
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depletion of 4 He is complete. Therefore, for r < 10, breakup no longer contributes to the 
cooling of the disk, its contribution being restricted to the viscosity through the neutrons 
already produced. 

To underline the pion contribution relative to the neutron production let us calculate it 
at the point where the ion temperature is largest, i.e., r — 1. Calculating physical variables 
at this point and inserting in eq. (26) yields 

y + 2 + 20/19 y + - 0.364 cT 2 Tf 2 t(w) , (38) 

and a is now (eq. (33)) 

a? = 730 ( 1 + 19 M T -1-85 _ 4 x i 3 T -3.7 (39) 

To see how a and y + behave in that region, we plot them in Figs. 6 and 7. 

In Fig. 6 we see that production of neutrons due to p (p 7r + n) p may be comparable to 
the contribution due to the breakup. 

Looking at eq. (35) we realize that the constraint on the accretion rate now reads 

Mi 7 < 0.183 Ti 1 - 85 . (40) 

As we have done previously, we now investigate the effects of cooling and pair production 
and, for that, we reduce our system of equations to only two equations involving T e and Tj. 
The thermal equilibrium and pair equilibrium equations turn to, respectively 

(20/19 + y + f - 3.86 a' 2 T~ 2 t{w) - 

1.1 a -1 T' 05 T e ~ 2 - 5 (1 + 0.41T e a5 ) (20/19 + y+) = 

and 

A 2 - 1 - 204 T 2 (1 + 0.41T e a5 ) 2 e" 11 - 87 ^ T e ~ 10 = , (41) 
where A is now given by 

A = 0.1975 a TjT e _1 (20/19 + y+) — y + , (42) 

and y + and a are given respectively by eqs. (38) and (39). 

In Fig. 8 we plot the solutions to the thermal and pair equilibrium equations. 

To see how the solutions of these equations behave as a function of the accretion rate, 
we have plotted them for M 17 = 0.9 in Fig. 9. There it can be seen the emergence of two 
solutions instead of only one as in the M 17 = 1.0 case. 
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In Fig. 10 we plot (1 + 2z) as a function of T; L . In contrast to Fig. 5, we have now two 
solutions, both with z < 1.0. 

Finally, to see the importance of nuclear cooling due to n + production, we plot q = 
Q pp /Q e i in Fig. 11. 

8. Analysis and conclusions 

We don't know by sure what physical processes may be operating in the outermost 
parts of the disk generating viscosity and switching accretion on. Putting this flaw aside, we 
realize that under certain conditions protons and electrons are out of thermal equilibrium, 
the ion temperature being much greater and close to the virial one and, for a matter of self- 
consistency, we are compelled to consider some nuclear reactions, since we have available 
energy up to twice the proton rest mass. However, besides energy considerations, we had 
a very strong motivation to consider 4 He breakup and p{pwK + )p. This is related to the 
viscosity problem itself: these reactions produce neutrons whose collisions with the accreting 
matter may be a source of the required viscosity to drive accretion. We have considered these 
reactions in an environment of protons, electrons, neutrons, photons and pairs. Concerning 
radiative cooling and pair production, we have made some simplifications that overestimated 
these processes and, as a consequence, underestimated the nuclear reaction rate because the 
total energy available is constant. The treatment we have adopted to calculate the nuclear 
reaction rate consists of an improvement over a previous one (Meirelles 1993) due to the 
abandon of the assumption of equality between nuclear reaction and dynamical times, as 
well as for the inclusion of both pion and pair production. To emphasize the role of the 
nuclear reactions, as far as the requirement of huge ion temperatures needed to ignite pion 
production, we have assumed a maximally synchronously rotating Kerr black hole, the inner 
radius of the disk equal to its horizon, R^. The region of the disk we were concerned is 
the one extending from the horizon to about 20 R h . In the outer parts of the inner region, 
the contribution of the reaction p(p, mr + ) p to the production of neutrons is negligible, this 
being dominated by 4 He breakup, which we use as a kind of boundary condition. Neutrons 
produced by this reaction contribute to the viscosity and the plasma heats up as we approach 
the hole, making possible the production of pions which, in turn, increase the viscosity by 
concomitant production of neutrons. 

The radiative cooling of the disk was assumed to be due to unsaturated inverse Comp- 
tonization of soft external photons impinging upon the disk. Since we have not considered 
the reaction p(p, pn°)p, we have not taken into account the radiative cooling due to photons 
coming from the n° decay. Part of those soft photons are upscattered in energy, reaching 
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the Wien peak where they interact with themselves producing pairs. Radiative cooling and 
pair production decrease the electron temperature, an effect that hinders the efficiency of 
nuclear reactions, even if the ion temperature is great enough for them to occur. When T e 
decreases, the number of interacting electrons with one proton increases sharply, due to less 
shielding. Nuclear forces only act at short distances and the high ion temperature needed 
to overcome the Coulomb barrier also hinders the occurrence of nuclear events due to very 
short time in the nuclear range. As a matter of fact, reaction rate due to electromagnetic 
interactions prevails over the rate due to nuclear interactions if T e satisfies 

T e < 28.56 T/ /3 . (43) 

However, the energy transferred in the (electromagnetic) scattering will be in the KeV range, 
while for the nuclear reaction it will be in the MeV range. Therefore, the previous inequality 
changes roughly to 

T e < 0.2856 T; 1/3 , (44) 

when one also considers the energy transferred in the event. 

Aiming at the application to Cygnus X-l and related sources, we have set My? = 1 and 
Myj = 0.9 which are much less than the critical value for Myj calculated close to the horizon, 
i.e., 

M{ 7 = 0.183 T; 185 . 

For both values of the accretion rate T e is larger than the value given by the ineq. (43). As a 
result, for every solution we have found, z < 1. This is surprising, since it is known that, in 
the absence of these reactions, one of the solutions for the two-temperature disk has z > 1. 
Another surprising result we have found concerns the sensitivity of the number of solutions 
to the value of the accretion rate. Calculating at r = 1 and r = 10, for M 17 = 1 we have 
found the solution is unique. Decreasing to M 17 = 0.9 there are two solutions. For these 
values of the accretion rate the neutron abundance in the inner region is fairly high, going 
from y n = 1/19 in the outer border, with total 4 He depletion, to about twice that value 
close to the horizon. This implies a quite large, i.e., large viscosity. This large viscosity, 
nevertheless, comes along with a high nuclear cooling, being comparable to the radiative 
cooling. 

It is known that the innermost region of the accretion disk is secular and thermally 
unstable (Shakura and Sunyaev 1976)in the absence of nuclear reactions. Taking into account 
these nuclear reactions, the disk behaves in the same manner at the very onset of these 
instabilities. However, it is worthwhile observing the behavior of the viscosity parameter 
with Tj, because two scenarios may emerge as far as the time evolution of the inner region is 
concerned, a starts growing, reaches a maximum somewhere between Tj = 15 and = 20, 
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then decreases and may become very small, even null, with increasing Tj. At that moment, 
accretion starts switching off and the cooling, that goes with a -2 , is practically instantaneous. 
Outside this region, however, some other mechanism for viscosity generation is still operating 
and accretion, there, keeps going on. As the matter reaches the border of the inner region 
it gets piled up there. As this piling up goes on, this region will be subject to several 
instabilities. Once one of these instabilities starts to grow it triggers accretion in the inner 
region. As the accretion proceeds, the accretion rate decreases till it reaches a value for which 
there are two steady states for the disk. Which one will the disk choose? At that moment, 
do both thermal and pair equilibrium hold in that region? A less drastic scenario is the one 
somehow similar to that proposed by Kusunose and Minishige (1991). Physical processes in 
the disk are characterized by time scales, the one chosen by the system being the least scale. 
As we have seen, thermal and pair equilibrium only hold under special circumstances. It 
may happen that due to the nuclear reactions and instabilities in the disk, the system may 
undergo a kind of limit cycle behavior around the upper solution. 

However, to have a better understanding of the time evolution of these systems we 
should make better treatment of both radiative cooling and pair production and a more 
detailed stability analysis of the problem as well. That is what we intend to do next, in a 
future contribution, taking also into account the reaction p {p pn°)p. 

C.L. Lima acknowledges support from FAPESP (project nr. 1998/065902) 
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Fig. 2. — The neutron abundance y n , as a function of the ion temperature, calculated at 
r = 10 for Mi 7 = 1.0. For Tj = 4.16 there is total 4 He depletion. 




4. — Relation T e — T { obtained through thermal equilibrium (lower curve) and through 
equilibrium (upper curve). 
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Fig. 5. — Relation between (1 + 2z) and Tj, for Mn = 1.0. The upper curve is obtained 
using Tj — T e data from the solution of the thermal equilibrium equation; the lower, using 
data from the pair equilibrium equation. 
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Fig. 7. — The viscosity parameter a taking pion production into account at r = 1. 
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Fig. 8. — The solution to the thermal equation (upper curve) and the pair equilibrium 
equation (lower curve). The curves are tangent at 7* 21. 
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Fig. 9. — The solutions to the thermal and pair equilibrium equations for M 17 = 0.9. 
Contrary to the previous figure, there are two solutions. 
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Fig. 10. — (1 + 2z) for M 17 = 0.9. The curve with larger values at the right is obtained 
using the data from the pair equation, the other using data from the thermal equation. 
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Fig. 11. — q = Qpp/Qei for M 17 = 0.9. The upper curve on the right is built with data from 
thermal equation, the other with data from the pair equation. For Tj < 12 the curve coming 
from the thermal data grows steeply (not shown here). 



